The Design of Random Surfaces with Specified Scattering Properties: 

Surfaces that Suppress Leakage 



A. A. Maradudin Q , I. Simonsen a - fc , T. A. Leskova c , E. R. Mendez d 
a Department of Physics and Astronomy and Institute for Surface and Interface Science, 
University of California, Irvine, CA 92697, U.S.A. 
b Department of Physics, The Norwegian University of Science and Technology, 

N-7491 Trondheim, Norway 
c Institute of Spectroscopy, Russian Academy of Sciences, Troitsk 142092, Russia 
Division de Fisica Aplicada, Centro de Investigacion Cientifica y de Educacion Superior de Ensenada, 
Apartado Postal 2732, Ensenada, Baja California 22800, Mexico 

We present a method for generating a one-dimensional random metal surface of finite length 
L that suppresses leakage, i.e. the roughness-induced conversion of a surface plasmon polariton 
propagating on it into volume electromagnetic waves in the vacuum above the surface. Perturbative 
and numerical simulation calculations carried out for surfaces generated in this way show that they 
indeed suppress leakage. 

I. INTRODUCTION 

The great majority of theoretical studies of scattering from randomly rough surfaces have been devoted to the 
direct problem, in which the surface profile function and its statistical properties are specified, and it is the angular 
distribution of the scattered intensity, and its polarization properties, that are sought. Comparatively little attention 
has been devoted to the inverse problem. In the usual formulation of this problem scattering data are provided by 
experimentalists, and the surface profile, or some statistical properties of it, are sought. This kind of inverse problem 
can be called a passive inverse problem, because the theorist has no control over the experimental data provided. 
In contrast, in the present work we consider a different type of inverse problem, one that we can characterize as an 
active inverse problem. In this type of problem, the angular distribution of the scattered intensity is specified, and 
the problem is to design, and ultimately to fabricate, a surface that generates that angular distribution. 

This active type of inverse problem has been studied even less than the passive type, yet it arises in a variety of 
contexts of both a technical and a basic physics nature. In this paper, we present an illustrative example of the design 
of a random surface with specified scattering properties that is prompted by the observation that as a surface plasmon 
polariton propagates across a randomly rough metal surface it continuously loses energy through its roughness-induced 
conversion into volume electromagnetic waves in the vacuum above the random surface, that propagate away from 
the surface. This leakage, as it is called, interferes with the determination of the Anderson localization length of 
the surface plasmon polariton by means of numerical simulation calculations, or experimental measurements, of its 
transmissivity as a function of the length of the random segment of the surface [1-3]. The use of a random surface in 
such studies that suppresses leakage therefore facilitates the investigation of the strong localization of surface plasmon 
polaritons by random surface roughness. 

In the approach to the suppression of leakage taken by Sornette and his colleagues [1,2], it was assumed that the 
random surface was not planar on average, but periodic, so that the dispersion curve of the surface plasmon polaritons 
supported by the mean surface displays a gap at the boundary of the one-dimensional first Brillouin zone defined by 
the period of the mean surface. Leakage should then either vanish or substantially decrease for the surface plasmon 
frequency at the band edge. However, this was not observed in the numerical simulation calculations of leakage carried 
^ ' out in Refs. 3 and 4. 

In this work we present an approach to designing a one-dimensional random surface that suppresses the leakage of a 
surface plasmon polariton as it propagates across it that differs from that proposed by Sornette et al. [1,2]. Although 
\ the power spectrum of the resulting surface is nonzero in a narrow range of wave numbers, that surface is not periodic 
on average. However, as with the surface proposed by Sornette et al., our surface is specific to the frequency of the 
surface plasmon polariton propagating across it: if that frequency is changed, a new surface has to be designed. 
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FIG. 1. The system studied in this paper. 



II. THE SCATTERED FIELD 



We study the scattering of a p-polarized surface plasmon polariton of frequency ut propagating in the xi -direction 
that is incident on a segment of a onc-dimcnsional randomly rough surface defined by the equation £3 = C( x i)- The 
surface profile function Ci x i) is assumed to be a single-valued function of x\ that is nonzero only in the interval 
— L/2 < xi < L/2 (Fig. 1). The region x 3 > ((xi) is vacuum; the region x 3 < ((xi) is a metal characterized 
by an isotropic, frequency-dependent, complex dielectric function e(ui) = ei(ui) + ie 2 (oj). We are interested in the 
frequency range in which e\(u)) < — l,e 2 (^) > 0, within which surface plasmon polaritons exist. We assume that the 
surface roughness is sufficiently weak that the surface profile function ((xi) satisfies the conditions for the validity 
of the Rayleigh hypothesis [5]. In this case the single nonzero component of the magnetic field in the vacuum region 
£3 > ((xi ) m ax can be written as the sum of the fields of the incident and scattered waves 
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In Eqs. (2.1)-(2.2) k{u) is the surface plasmon polariton wave number, 
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characterize the exponential decay of the field of the surface plasmon polariton with increasing distance from the 
interface into the vacuum and the metal, respectively. The functions R > (q, ui) and R < (q,to) are the scattering 
amplitudes of the surface plasmon polariton in the vacuum and in the metal, respectively, and 

a (q,u)) = (^j - q 2 ^j Reao(q,w) > 0,Jma o (?,w) > (2.6a) 

a(q, uj) = f e(uj)— - q 2 J Rea(q, oj) > 0, Ima(q, oj) > 0. (2.6b) 

The scattering amplitude R > (q, oj) satisfies the reduced Rayleigh equation [6] 
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This equation will be solved perturbatively and numerically. 

III. THE TOTAL SCATTERED POWER 

From Eq. (2.1) we see that the scattered field in the vacuum region can be written in the form 

H>( Xi ,x s \oj) sc = / p-R>(q,u)e^ +iao ^ x '. (3.1) 
The total scattered power, normalized by the total power in the incident surface plasmon polariton, is 

S(oj) = (3.2) 

where 
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with L 2 the length of the surface along the 22-axis. The scattering angle 9 S , measured clockwise from the x 3 -axis, is 
related to the wavenumber q by q = (oj/c) sin 9 S . We therefore find that 

S(oj) = 1 ^- £{UJ))3/2 r d6 sC os 2 9 s R> (-sm0„u) \ (3.4) 

Since the integrand in Eq. (3.3a) (and in Eq. (3.4)) is non-negative, we see that the only way in which leakage can 
be suppressed, i.e. the only way in which P sc can be made to vanish, is to design a one-dimensional random surface 
for which the amplitude R > (q,oj) is identically zero for — (uj/c) < q < {oj/c). In the next section we introduce one 
such surface. 
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IV. THE RANDOM SURFACE 



We write the surface profile function C(x\) in the form 

C(x 1 )=T(x 1 )s(x 1 ), 



(4.1) 



where s(x\) is a single-valued function of x\ that is diffcrentiable and constitutes a stationary, zero-mean, Gaussian 
random process defined by 
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where the angle brackets denote an average over the ensemble of realizations of s(xi). The function T(xi) serves to 
restrict the nonzero values of s(xi) to the interval — L/2 < x\ < L/2. One form T{x\) can have is 



T(x 1 ) = 0(^+x 1 )0(^-x 1 ), 

where 6{x\) is the Heaviside unit step function. A smoother, diffcrentiable version of T(x\) is provided by 

1 + cosh \f3L 



T(xi) = 
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where the parameter controls the range of x\ values over which T(x\) decreases from 1 to 0. A value of /3 given by 
100/L cuts off s(xi) smoothly. 

The power spectrum of s(xi), g(\Q\), is defined by 

(4.5) 
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A surface that suppresses leakage is defined by the power spectrum (Fig. 2) 
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FIG. 2. The power spectrum of the surface roughness assumed in this work 



That a surface characterized by the power spectrum ([hfj) suppresses leakage can be seen from the following argu- 
ment. The incident surface plasmon polariton has a wave number whose real part is ki(u>). After its first interaction 
with the surface roughness the real part of its wave number will lie in the two intervals (3ki(u>) — Afc, 3k\(u)) + Afc) 
and (— ki(u>) — Afc, — ki{ui) + Afc). This is because the wave numbers in the spectrum of the surface roughness with 
which kx(u>) can combine lie in the intervals (2k%(u)) — Ak,2kx(w) + Afc) and (— 2ki(u>) — Afc, — 2k\(uj) + Afc). For 
the same reason, after its second interaction with the surface roughness the real part of the wave number of the sur- 
face plasmon polariton will lie in the three intervals (5fci(w) — 2 Afc, 5fci(w) + 2Afc), (ki(u>) — 2Afc, fci(w) + 2Afc), and 
(— 3fci(u>) — 2Afc, — 3fci(w) + 2Afc). After three interactions with the surface roughness the real part of its wave number 
will lie in the four intervals (7fci(w)-3Afc, 7fci(w)+3Afc), (3fci(w)-3Afc, 3fci(w)+3Afc), (-fci(w)-3Afc, -fci (w)+3Afc), 
and (— 5fci(cj) — 3Afc, — 5fci(w) + 3Afc), and so on. Thus, for example, if — fci(cj) + 3Afc < — (cj/c), so that 
Afc < ^(ki(uj)~- (ui/c)), after three scattering processes the surface plasmon polariton will not have been converted into 
volume electromagnetic waves. In general, if we wish the surface plasmon polariton to scatter n times from the surface 
roughness without being converted into volume electromagnetic waves, we must require that Afc < i(fci(u;) — (uj/c)). 
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FIG. 3. A perturbative result for (uj/c) 2 (\R> (q,uj)\ 2 } as a function of (cq/ui) for a surface plasmon polariton, whose frequency 
corresponds to a vacuum wavelength A = 457.9 nm, propagating across a random silver surface (e(w) = —7.5 + i0.24) of length 
L — 20A, defined by the power spectrum (4.6) with Afc = 0.3(fci(w) — (u/c)), 8 = 10 nm, and Eq. (4.3). 
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FIG. 4. A computer simulation result for (lo/c) 2 (\R > (q, uj) | 2 ) as a function of (cq/ui), plotted on (a) a linear-linear scale and 
(b) a linear-log scale, for a surface plasmon polariton, whose frequency corresponds to a vacuum wavelength A = 457.9 nm, 
propagating across a random silver surface (e(co) — — 7.5 + i0.24) of length L = 20A, defined by the power spectrum (4.6) with 
Afc — 0.3(fci(cj) — (u/c)), 8 — 30 nm, and Eq. (4.4). The results for 50 realizations of the surface profile function were averaged 
to obtain the results plotted in this figure. 

To show that this approach to suppressing leakage works, in Fig. 3 we present a plot of (u/c) 2 (\R > (q, u)\ 2 ) as a 
function of cq/u for a silver surface characterized by the power spectrum (4.6) with Afc — 0.3(ki(u) — (u/c)) and 
S = 10 nm. The frequency of the surface plasmon polariton corresponds to a vacuum wavelength of A = 457.9 nm, 
and the dielectric function of silver at this frequency is e(u) — —7.5 + i0.24. The calculation of (\R > (q,u)\ 2 ) was 
carried out as an expansion in powers of the surface profile function through terms of fourth order. For this i? > (<;, u) 
was calculated as an expansion in powers of C( x i) through terms of third order (there is no zero order term). The 
expression (4.3) was used for T(xi), and the length of the random segment of the metal surface was L = 20A. The 
peak in (u/c) 2 (\R > (q,u)\ 2 ) for q < —(ui/c) arises from the terms in the expansion of R > (q,u) of first and third 
orders in C( x i); that for q > (u/c) arises from the second-order term. We see from this figure that (\R > (q, u)\ 2 ) 
indeed vanishes for — (u/c) < q < (u/c). With our choice of Afc = 0.3(fci(w) — (u/c)), a calculation of R y (q,u) to, 
say, fourth order in ((xi) would have yielded a small nonzero contribution to (\R > (q,u)\ 2 ) in the radiative region 
— (u/c) < q < (u/c), but even that could be suppressed by choosing Afc < 0.25(fci(w) — (u/c)). 

In Fig. 4 we show that the same result is obtained if we increase the roughness of the surface by increasing 5 from 
10 nm to 30 nm, keeping the same values for the remaining material and roughness parameters that were used in 
obtaining the results plotted in Fig. 3. In this case a computer simulation approach based on the numerical solution 
of Eq. (2.7) was used. The expression (4.4) was used for T(xi), and the results for 50 realizations of the surface 
profile function were averaged to obtain the results plotted in this figure, in Fig. 4b four peaks are easily seen. 
They correspond to the real parts of the wavenumbers of the scattered surface plasmon polaritons resulting from the 
scattering of an incident surface plasmon polariton of wavenumber k(u) = k\(u) + ik 2 (u) = (1.0741 + «0.0026)w/c. 
Earlier in this section we predicted that the real parts of these wavenumbers should be in the vicinity of ±fci(w), 
±3k\(u), etc. These predictions fit very well with what can be observed from Fig. 4b, even though the peaks at 
±5k\(u), and higher, cannot be distinguished from the noisy background due to the low number of samples on which 
this figure is based. It should be observed that the peak at k\(u) is stronger then the one at 3k\(u), even though 
the latter is of first order (in perturbation theory) while the former is of second order. We believe that this is related 
to the fact that both the surface plasmon polaritons resulting from a single interaction with the rough surface can 
be scattered into surface polaritons of (real) wavenumber k\(u), which in this sense gives it a multiplicity of two. 
The surface polaritons at 3fci(w) can, in addition to the single scattering contribution, only be reached by scattering 
processes of order three or higher, which are all weak for the roughness assumed here. 

Thus, the method presented in this work provides an effective way of producing random surfaces that suppress 
leakage. 
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V. CONCLUSIONS 



In this paper we have presented an approach to generating a one-dimensional random surface that suppresses 
leakage. We have shown by perturbative and numerical simulation calculations that surfaces generated in this way 
indeed possess this property. 
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